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A simple calculation method for photon statistics of frequency-filtered fields is proposed. This 
method, based on eigenvalue decompositions of superoperators, allows us to study effects on the 
photon statistics of spectral filtering by various types of filters, such as Gaussian and rectangular 
filters as well as Lorentzian filters, which is not possible by conventional approaches. As an example, 
this method is applied to a simulation of quantum dot single-photon emitters, where we found the 
efficient choice of the filter types to have pure single photons depends on the excitation conditions 
i.e. incoherent or coherent (and resonant) excitations. 
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I. INTRODUCTION 


High quality single photon (SP) sonrces, which emit 
one photon at a time with high parity and high rate, are 
essential for the realistic and reliable application to qnan- 
tum information science and technologies Qi- As effi¬ 
cient solid-state SP sources, semiconductor quantum dots 
(QDs) are promising candidate systems in solids and have 
been attracting attention for a nnmber of advantages: the 
well-defined quantized states, high controllability in the 
emission wavelength, high brightness even enhanced by 
embedding them in nanocavities, the emission-site con¬ 
trollability, and possible cnrrent injection operations 

El- 

However, in QD SP emitters, a number of emission 
lines are typically present due to the multiple transition 
levels and also to other QDs in a sample, degrading the 
SP purity. To avoid the degradation, spectral filters (Fig. 
[T]) are usually used for selecting the relevant emission e.g. 
an exciton emission, and filter ont the spectrally sepa¬ 
rated irrelevant emissions, snch as, the biexciton-exciton 
emission 0 and charged excitons. The role of nsing 
a spectral filter is to prevent the detection of irrelevant 
emissions spectrally separated in the frequency domain. 
At the same time, dne to the frequency-time nncertainty, 
using narrow spectral filters inevitably widens the detec¬ 
tion field in the time domain, leading to degradation in 
the time resolntion and also in the pnrity of the SP emis¬ 
sions. In this way, spectral filter modifies the filtered held 
both in the freqnency and time domains, and therefore, 
the hltering effect on the photon statistics (which basi¬ 
cally is a mnltiple-time correlation function given in the 
time domain) is not so simply understood, especially for 
quantum emitters. 

Theoretical and experimental stndies of hltering effects 
on the photon statistics have recently been attracting 
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attention [l5l - 0 . These stndies were triggered by the 
development in the theoretical treatment, a versatile cal- 
cnlation method proposed by E. del Valle et al. 0- In 
this method, the spectral hltering process is effectively 
replaced with the inclusion of probe systems coupling 
weakly to the system. A g reat advantage of this method 
over the former theory |20l422l | is that the complication in 
calculation coming from the time orderings of operators 
can be avoided, allowing the calculation of higher-order 
correlation functions (n > 3). However, in this method, 
the type of spectral hlters the method can treat is re¬ 
stricted to Lorentzian hlters, since the spectral hlter is 
mimicked by a Lorentzian density of states of a probe 
system under the Markov decay process. The former an¬ 
alytic, but approximate, approach [2l| also treats only 
Lorentzian hlters due to its simplicity. Therefore, the ef¬ 
fect of spectral hltering on the photon statistics has been 
investigated only for Lorentzian hlters so far. 

In this paper, we propose a simple calculation method 
that allows for theoretical treatment of a variety of spec¬ 
tral hlters in order to deepen the understanding of hlter 
effect on the photon statistics. In Sec. im we introduce 
the calculation method based on superoperator eigen- 


quantum emitter 

filter 


HBT set up 



FIG. 1: Quantum emitter and detection system. The emission 
dynamics of the quantum emitter is given by the quantum 
master equation, ■§ip = Cp. The effect of spectral filter is 
described by the correlation function fir), and the photons 
passed through the Hlter enter into the HBT setup [T^ for 
the correlation measurement. 
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value decomposition, and find the exact expressions for 
the second-order correlation function for Gaussian and 
rectangular filters as well as the Lorentzian filters. As 
far as we know, this is the first time to apply the su¬ 
peroperator eigenvalue decomposition technique to the 
higher-order correlation functions, whereas it has been 
applied to the calculation of the first-order correlation 
functions e.g. in a calculation of the Mdssbauer spec¬ 
tra [ 2 ^. While our method allows for the treatment of 
types of filters, it directly treats the operator ordering 
problem, and thus the difficulty in calculating high order 
correlation functions is not removed. In this sense, our 
method is complementary to the previous theory [l5l| . In 
Sec. mil as an example, we show a numerical simulation 
applied to QD SP emitter systems, where we found the 
efficient choice of the filter types for purifying the single 
photons depends on the excitation conditions i.e. inco¬ 
herent or coherent (and resonant) excitations. 

We note that the effect of the background noise which 
is not related to the system dynamics is out of the scope 
in the theoretical framework. Throughout the paper, we 
set h = 1 for simplicity unless otherwise specified. 


II. SUPEROPERATOR EIGENVALUE 
DECOMPOSITION METHOD FOR PHOTON 
STATISTICS OF FREQUENCY FILTERED 
FIELDS 

A. Definition of the problem 

Here, we will define the problems to solve. The sys¬ 
tem we consider consists of an emitter, a spectral filter, 
and detection system as shown in Fig. [T] Photons emit¬ 
ted from the quantum emitter are detected by Hanbury- 
Brown Twiss (HBT) setups [l^ for the second-order in¬ 
tensity correlation measurement, and before the detec¬ 
tion, the photons passed through a spectral filter whose 
response is described by the filter correlation function 
/(r). Alternatively, the detection can be performed with 
high-speed streak camera with high time resolution (less 
than a few picoseconds), which is now becoming a pow¬ 
erful detection system for the study of photo-counting 
statistics [ 2 ^. Throughout the paper, we assume the ef¬ 
fect of back reflection at the filter surface on the emitter 
system can be neglected [ 2 ^. In this case, the dynamics 
of the quantum emitter and the emission field is given by 
the Hamiltonian, for the emitter, Lindblad-type su¬ 
peroperators, Crj, for decay and pump processes labeled 
by r] with the rates, 7 ,^, and the resulting quantum master 
equation [ 2 ^, 

^^p = i[p,H]+^')r^Cr^p = Cp. (1) 

77 

The emission field operator, at the exact emission 
time, t, is given by the Heisenberg operator (account¬ 
ing for the system dynamics except for the filter and de- 



FIG. 2: Schematics of filter function: (a) /(r) in time domain 
and (b) in frequency domain. 


tection systems), and the frequency-filtered field 

Ep{t) to be detected at a time t is given by 

POO 

EFit)= f{T)E~{t-T)dT, 

Jo 

pOO 

Et{t)= r(r)A+(t-r)dr, (2) 

Jo 

The time region of the integration is physically restricted 
to T > 0 by the causality, and the correlation function 
in the time domain, /(r), has the peak at a delay time, 
T = Td, corresponding to the optical path length between 
the emitter and the filter, and it has a width Tc corre¬ 
sponding to the hlter correlation time (Fig. [2] (a)). The 
filter function F{ui) in the frequency domain, is centered 
lof and has a bandwidth A roughly equal to the in¬ 
verse of the correlation time, A ^ 1 /tc (Fig. (b)). The 
equations 0 and m are the most general expression for 
system dynamics and the filtered emission field, and thus, 
can be directly applied to any emitters and any filters. 
For example, in case of a resonantly driven two-level atom 
(transition energy wa, Rabi frequency Hr, and the laser 
frequency Wi), H = 

the spontaneous emission decay is included with the rate 
dv ~ and superoperator in standard nota¬ 

tion, and the emission field is E^{t) = 

The n-th order normalized intensity correlation func¬ 
tion to be evaluated is then given by 

,tn) 

_ {T+T-Epiti) ■ ■ ■ EF{tn)Ep{tn) • ■ ■ Ep{ti)) 

(3) 

where 7k and 7+ are time ordering and anti-ordering 
superoperators working on the Heisenberg annihilation 
and creation operators, respectively. The brackets mean 
the ensemble statistical average over the emitter states, 
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and mathematically given by taking the trace after mul¬ 
tiplying by the density matrix of the emitter system, 
(O) = Tr(0/5). In addition, from Eq. ([2]), it is necessary 
in calculating Eq. ([3]) to evaluate the operator products 
with different time arguments by using quantum regres¬ 
sion theorem [2^. The aim of this paper is to give a 
simple calculation method for the correlation function in 


Eq. (I3|). 

The time ordering operation has to be taken into ac¬ 
count in Eq. ([3|) when the effect of back reflection by the 
filter is negligible as we assumed here [2^. However, the 
operation makes the calculation of the n-th order cor¬ 
relation function of large n (like n > 3) complicated. 
For = • • • = t„), the number of time arguments 

{tj with j = 1 , • • • n for fields and j = n -I- 1, • • • 2n 
for E~ fields) are 2n. The number of different time or¬ 
derings in the integration is reduced to (2n)!/(n!)^ by a 
symmetry argument; The number is reduced from (2n)! 
to (2n)!/(n!)^ since the time ordering operators 7+ and 
71 sort the product of E+ fields and E~ fields of nl dif¬ 
ferent orderings, respectively, into one exclusive ordering. 
Therefore, the number of terms with different time order¬ 
ings for the n-th order correlation function amounts to 
six for n = 2, twenty for n = 3, and seventy for n = 4 211. 

In our approach given below, we will finally obtain 
the analytic expression for the correlation function in 
Eq. ([3l), whereas the time ordering process is directly 
treated, hence the difficulty is not removed. Therefore, 
the higher order correlation function with n > 4 is too 
computationally expensive. In this sense, this method 
is limited to the application to the correlation functions 
with n < 3 in realistic calculation, while the recently 
proposed method [HI can avoid the complicated time or¬ 
dering operation to be able to simulate photon statistics 
to the higher order. 

However, as mentioned in the introduction, this 
method allows us to have analytic results for general 
types of the filter function /(r), whereas the previous 
method [l5| can treat only the Lorentzian filters. In this 
sense, our method is complementary to the other meth¬ 
ods ill, [HI, and this makes possible the comparison of 
the efficiency in optimizing the photon statistics for dif¬ 
ferent types of filters. 


B. Superoperator eigenvalue decomposition 
method 


written in a linear equation -^O = LO after re¬ 
forming the operator O into a vector form O = 

- , Ol.AT,., 02,1, • • • ,02,N^,-“ jOat,,,!,--- ,ONc,Nc) 
The length of O is N^, and the Liouvillian matrix L 
(originally the superoperaor £) has dimension, x N^. 
Therefore, L has eigenvalues, H, which are in general 
complex values with Re(H) < 0. The corresponding 
right and left eigenvectors, vq and Uq, are defined here 
as 


Lvq = flvQ, MqL = UqQ. (4) 

Therefore, if the eigenvalues are non-degenerate, the op¬ 
erator in vector form, O, is decomposed into the eigen¬ 
vectors: 


d = ^0(0)un, (5) 

Q 


where 0(12) = (mq • O )/ (mq -vq). From the above expres¬ 
sion, we obtain the spectral decomposition of the vector 
by d = Hn[^](^) [d](H) = C(0)un, whose matrix 

form is written as O = Hn[d](f^)- The merit of using 
the eigenvalue decomposition is that the time evolution 
of the operators are explicitly given by 

d{t)^e^*[6]=J2[dme^K (6) 

a 


Average filtered-field intensity — Now, we apply the 
eigenvalue decomposition technique to the filtered-field 

intensity, {E+{tj)E~{tj)) = Tr (^E+{tj)E~{tj)pss^ in 

Eq. Q, whereas the system is assumed to be in the 
steady state p{tj) = pss- Inserting Eq. Q, we have 


fOO poo 


{Ep{tj)EF{tj)) — 


/o ^0 


dridr2 /*(n)/('r2) 
xTr (^E+{tj - Ti)E~{tj - T 2 )pss'^ 

= [[ dridr2 r(n)/(T-2)Tr(d+e^"-[d-/5ss]) 

+ [[ dTidT2 /*(ri)/(r2)Tr (£'“e^^i"[/3ss£'+]') , 


(7) 


with Tij = Ti — Tj. Applying Eq. ([51), the filtered-field 
intensity is expressed in the form 


Here, we will introduce our method based on super¬ 
operator eigenvalue decompositions. As an example for 
n = 2, we will obtain a general expression for g^^^(T = 
ti — t2 = 0). 

Superoperator eigenvalue decomposition — According 
to the quantum regression formula [2^, different¬ 
time correlation functions can be calculated by the 
same equation as the density matrix equation in 
Eq. (P). Thus, any time dependent operator 0{t) 
satisfies -iO = £0. The matrix equation can be 


{E+{t,)Ep{t,)) = Y,smg{n) + s*{n)q*{n), (s) 

n 

where we noticed that the first and second terms in the 
r.h.s. of Eq. ([7]) are the conjugate pairs. The coefficients 
are given by 

s(H) = [[ dridr 2 /*(Ti)/(T 2 )e^^=\ (9) 

T2 >Ti >0 

g(H) = Tr(d+[d-pss](H)) . (10) 
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This is the general form of the superoperator eigenvalue 
decomposition for the filtered-field intensity. If the filter 
bandwidth A is set small, the intensity 
as a function of the central frequency wj’ is the emission 
spectrum. This eigenvalue decomposition method was 
previously applied to the calculation of the Mossbauer 
spectra as a first-order correlation function. The 
method shown here is essentially the same as that shown 
in the paper. However, we will now apply this method 
to the second-order correlation function. 

Average filtered-field intensity correlation — Next, we 
compute the second-order correlation function at zero de¬ 
lay assuming the steady state. 


{r+r-E+{ti)E+{t 2 )Efi{t 2 )Efi{ti)) 


X 


r4>r3>0 

Tr (^E+{-Ti)E+{-T2)E~{-T3)E~{-T4)pss^(fi^) 


where dr"^ = dTidr 2 dr 3 dr 4 . The region of the four-fold 
integration is divided into six regions with different time 
orderings: (i) r 2 < T 3 < r 4 < n, (ii) T 2 < r 3 < n < T 4 , 
(iii) r 2 < Ti < r 3 < r 4 , (iv) T 3 < r 2 < n < r 4 , (v) 
T 3 < T 2 < T 4 < ti, (vi) T 3 < T 4 < T 2 < Ti. Since the 
contributions from (i) and (iv), (ii) and (v), and (iii) and 
(vi), are complex conjugate pairs respectively, we have 
only to compute the integration over (i), (ii), and (iii). 

The second-order correlation function in Eq. (EH), 
At:{e+{-ti)E+{-T 2 )E~{-T 3 )E~{-T 4 )pss^, is ex¬ 
pressed, by using quantum regression theorem and the 
superoperator eigenvlue decomposition, as 


= TviE+e^'^^^ 






PssE^ 


= ^ exp(HiTi 4 -I- ^ 2 x 43 -I- H 3 T 32 ) (12) 

r2i ,^2 5^3 


xTr 




E- 


E- 


PssE^ (f^i) (^^ 2 ) (^^ 3 )) 


for (i) T 2 < T 3 < T 4 < Ti. Similarly, it is 


= ^ exp(Hir4i-I-r22ri3-I-H3r32) (13) 


cTr(£;+ 


E- 


E pss 


(Hi).E+] (H 2 )] (SI 3 )) 


for (ii) T 2 < T 3 < Ti < Ti, and 


= ^ exp(Hir43 -I- H 2 T 31 -I- H 3 T 12 ) (14) 

r2i ,r22 )^^3 


xTr 


[e- 


E- 


E PSS 


(Hi)] {n 2 )E+] (H 3 )) 


for (iii) T 2 < ti < T 3 < T 4 . By inserting the Eqs. m- 
(fTTl) into Eq. EH), we obtain a general expression 


(r+r_ E+ ih)E+ (t2)Efi it2)Efi (ti)) 

ti—t2 

iii 

= 2Re E E Zfe(Hi,H 2 ,H 3 ) 0 fe(Hi,H 2 ,H 3 ), (15) 

k—i fii,02)^3 


where 


Zi = 2^11dr4r(Tl)r(T2)/(T3)/(T4) 

xexp (H 1 T 14 -I-H 2 r 43 -I-H 3 T 32 ), (16) 

= 2 ' //drV*(ri)/*(r2)/(r3)/(r4) 

xexp(HiT4i-|-H2ri3-|-H3T32), (17) 

Ziii = 2^ Iff[ clrV*(Ti)/*(r2)/(r3)/(r4) 

J J J J (iii) 

xexp(HiT43 + H 2 T 31 -I-H3ri2), (18) 


and 


0; = Tr {e+ 
0ii = Tr {e+ 
0iii = Tr (e+ 


E- 


E- 


pssE+] (Hi)] (H 2 )] (H 3 )) , (19) 


E- 

E- 


E pss 
E~pss 


(Hi) {n 2 )E+ (H 3 ) , (20) 


(Hi)] (H 2 )E+] (H 3 )) . (21) 


With the general decomposed expression, the effect of the 
spectral filtering on the second-order correlation function 
enters only through Zfc(Hi, H 2 , H 3 ) and s(H). Therefore, 
they can be regarded as response functions of the system 
in which the filter response is convolved. 

We should mention here the case of short correlation 
time Tc filters (tc is defined through /(r Tc) = 0 ), 
which should correspond to an unfiltered case. If Tc is 
much shorter than the time scale of system dynamics, we 
can put exp(HiTj) = 1 for Zk in Eqs. (fT6l) - (IT8l) and for s 
in Eq. dH]). In this case (tc —>■ 0), s and ior k = i — iii 
are independent on H, Hi,H 2 , and H 3 . Then, bj using 
= {E+E+E-E-), = {E+E-), 

we safely find that the expression for the normalized cor¬ 
relation function is reduced to be that of the unfiltered 
field, 




{E+E+E-E-) 

{E+E-)^ 


( 22 ) 


s and Zk for Lorentzian/Gaussian/rectangular 
filters — As the typical examples, the above general 
expression is applied to three types of filters, Lorentzian, 
Gaussian, and rectangular filters, to obtain the explicit 
analytic forms for s and Zk here. In the calculation, 
we assume for simplicity that the time delay of the 
filter response, in Fig. [2] (a), is much larger than 
the correlation time, Tc, and in addition, the system 
is assumed to be in the steady state. Under this 
assumption, we will change the time variables from r 
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to T + Td and approximately change the lower limit of 
the time integration from 0 to —~ —oo. With this 
change, the range of the integration for s(n) is replaced 
by —oo < Ti < T 2 < oo in Eq. Similarly, for Zk 
in Eqs. (fT 6 ]) - (fT 8 l) . the time range of the integration is 
replaced by — oo < Tj < oo while the ordering among ti, 
T 2 , T 3 , and T 4 are unchanged. 

The Lorentzian filter is the simplest example to per¬ 
form the time integration to give s (= s^) and 
(= Z^), since the correlation function of Lorentzian filter 
/(r) = fL{T) is an exponential [ 2 l|, 

fLir) = A 6 »(r)exp ((-A - ta;F)r), (23) 


where 9{x) is the Heaviside step function. Inserting this 
and after straightforward integrations, we find that they 
are given by simple polynomial fractions. 


s^in) = - 

= 


A/2 


iio -j- A — H 
A A 


A 


A — iujF ~ Hi 2A — H 2 3A -I- itop ~ H 3 ’ 


zk = 


zh = 


A 


A 


A 


A -l- iojF — Hi 2A — H 2 3A -I- iiop — H 3 ’ 


A 


A 


A 


(24) 

(25) 

(26) 


A + iujp — 2A + i2L0F — 1^2 3A H- icvp — ^^3 

(27) 


For Lorentzian filters, the time integration for correlation 
functions gives the products of the transfer function, and 
therefore, analytic time integration upto the arbitrarily 
high orders is possible. Due to the simplicity, the pho¬ 
ton statistics of filtered fields has been studied only for 
Lorentzian filters [il,[2l|. However, as shown in the next 
section, when the time scale of the system dynamics is 
comparable to r^, which we sometimes face in state-of- 
the-art quantum emitters, the best choice of the filter 
type is essential. Therefore, the photon statistics of the 
field filtered by other types of filters should be necessary. 
Here we just show the results for Gaussian and rectan¬ 
gular filters (but the method can be applied to arbitrary 
hlter function). 

For Gaussian filters, the correlation functions, /(t) 
(= /g( 7 ‘)) in the time domain and F{uj) (= Fg{uj)) in 
frequency domain are, 

fair) = ^exp (-(At )2 - , (28) 



where the Fourier transform is defined by F{u}) = 
(27r)“^ f /(T)exp(ia;T)dr. For this filter, s(H) (= s‘^(H)) 
is given by 

■s‘^(H) = iexp (?/(H)2) (1-herf (?/(H))), (30) 


TABLE I: Coefficients for of Gaussian filters, Eq. (I31II . 


k 

Ak 

Bk 

Ck 

i 

Hi — 

— 2iujp — +f 23 

— 

2A 

2A 

2A 

ii 

— liuj p \ —^^2 ”l"^3 

— Qi +^2.3 

— 

2A 

2A 

2A 

iii 

Hi — 

— +^3 

2iu)p —0-2 

2A 

2A 

2A 


TABLE II: Coefficients for ZJ of rectangular filters, Eq. (I35II . 


k _Ofc_dfc_ Jk _ 

i ^^+i0 u^F+\-in^ 

-u:p-in3 -q -u:p+X-in^ 


where y(H) = (H — iojF)/{V2\) and erf(a;) is the Gauss 
error function. For the second-order correlation function, 
we obtained an analytic expression for Z^ (= Z^), 


ZF - —e^k+Bl+cl 


pOO 

/ (I-erf(z-fiC'fc)) 

Jo 


X (erf(z -I- Bk) - erf(-z -|- Bk)) dz. 


(31) 


whose coefficients, Ak, Bk, and Ck, are given in Table HI 
For rectangular filters, filter correlation functions, /(r) 
(= frij)) and F{ijj) (= Fr{uj)) are given by 


frir) 

. sin(Ar) 

= exp[—iujFr) -, 

(32) 


TTT 


Fr{w) 

= ^^(^ -\uj - wf )■ 

(33) 


For this filter, s(H) (= s’'(H)) is found to be 


s’'(H) 


I j /wf + A-|-iH — i0\ 
27ri ywF — A-|-fH — iOy ’ 


(34) 


where the infinitesimally small positive number, 0 , is in¬ 
troduced for the analytic continuation of the logarithmic 
function (which is essential in case Re(H) = 0). The an¬ 
alytic expression for Zk (= Z^) is also found by inserting 
Eq. (|32|l into Eqs. (fT 6 )l - (|T^ and performing the integra¬ 
tion. 


i ((F(«F, FF ; 2) + </(aF:/3F;-2)) 

x(ln (2 - 7 fe) - ln(- 7 fe)) 

-$(a+, ^+, 7G 2) + $(afc , /3fc-, 7fc - 2; -2)^ , (35) 

where ± I, /3^ = /3fc ± 2, and the coefficients ak, 

Pk, and 7 fc are given in Table H. The functions, (j) and 4>, 
are defined with an analytically continued function of the 
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n-th order polylogarithm, Li„( 2 :) = X]m=i by 


(j){a, b;z) = — ln(z — b) ln(—a) 
'z-b\ 

+ Ll2 
dz. 


+ lrL(z — a) In 

$(a, 6,c;z) = / 
JO 


a — b 
(t>{a,b;z) 
0 


z — c 



(36) 

(37) 


In the evaluation of Z'[, we have carefully performed the 
multiple complex integrations since the contours cross 
branch cuts of the logarithmic functions. 

To summarize this section, a simple calculation 
method for photon statistics of filtered field, based on 
superoperator eigenvalue decomposition technique, was 
proposed and analytic expressions for s and Zk are ob¬ 
tained for the three types of filters : Lorentzian, Gaus¬ 
sian, and rectangular filters (as typical examples). For 
other types of filters, it will also be possible to find 
analytic expressions, although we will not go into fur¬ 
ther details here. Validity of our method is confirmed 
numerically by the perfect agreements with the other 
method for the case of Lorentzian filters, as shown 
in Sec. HID 


III. APPLICATION TO QUANTUM DOT 
SINGLE PHOTON EMITTERS 

Here, we take QDs as an example of efficient SP emit¬ 
ters and apply the proposed eigenvalue decomposition 
method to a simulation of the photon statistics of the 
emission field filtered by Lorentzian, Gaussian, and rect¬ 
angular filters. 


A. QD SP emitters under incoherent pumping 

The model of the QD emitter system is the same as 
that used in our previous paper [Ij, . We consider 

the QD emitter states consisting of the electron-hole car¬ 
riers as shown in Fig. [3] (a). Among sixteen carrier con¬ 
figurations occupying the lowest energy levels, six charge- 
neutral configurations are taken into account: an empty 
state, |G), two bright exciton (BX) states, \BX1) = 
e|h||G) and \BX2) = e|h||G), two dark exciton (DX) 
states, \DX1) = e|/ij|G) and \DX2) = e|h||G), and 
a biexciton state, \XX) = e|e|h|/iJ|G), where ea and 

her (ej. and hj.) are annihilation (creation) operators of 
electrons and holes with spin a =t, 4 . in their respective 
lowest energy levels of the QD. The Hamiltonian of the 
QD emitter is 

H = u:xNtot-x\XX){XXl (38) 

where Ntot = X)<T=t number of 

excitons, x(= ujx — ojxx) is the biexciton binding en¬ 
ergy, and the fine structure splitting between the exciton 


(a)a = neutral, □ = charged 


G 

4 - 


— 

_ 

— 

— 


— <»■ 

BXl 

BX2 

DXl 

DX2 


— 

4 #- 


— 

4m 




—0- 

XX 

— «► 

4m 




(b) 


XX 


\ V 


BXl , BX2 DXl DX2 


FIG. 3: A neutral QD model [3, 1^ . 1^ . (a) Among 16 
electronic states at the QD ground levels, 6 neutral states 
with upto two excitons (G, BXl, BX2, DXl, DX2, and XX) 
are taken into account, (b) Incoherent pump (P) and decay 
processes (y^p, yl, 7s)- 


states is neglected. The following incoherent decay pro¬ 
cesses as are considered shown in Fig.[3](b): the decay of 
the injected electron-hole pairs is dominated by the spon¬ 
taneous emission (the rate ysp) [l^, the excitons suffer 
dephasing (with the rate Fph), and the spin flip of elec¬ 
trons and holes (with the rates y| and y^) results in the 
transitions between dark and bright exciton states with 
a rate ys (= 7 s + 7 s)- 

In this neutral QD, it was shown that XX emission at 
Lo = oJx ~ X is enhanced by incoherent XX excitation 
via DX states, and can strongly degrade the purity of SP 
emissions, especially in the case when the spin flip process 
is slower than the spontaneous emission (ys < y^p) 
Therefore, if the exciton emission at w = ojx is applied 
to a SP source, the XX emission must be effectively cut 
by using a spectral filter. Here, for the calculation of the 
emission properties, we define the emission field operator 
for the BX recombination by E~ = ^ah-a = 
and the central frequency of the filter wp is set luf = oJx- 

In Fig. |4] ( a), the emission spectrum is shown for a 
situation (x = 2 meV, Fph = 20 /reV, l/ys = 10 ns, 
l/ 7 sp=l ns). In the same figure, the filter functions in the 
frequency domain \F{uj)\'^ are also shown for Lorentzian, 
Gaussian, and rectangular filters with the bandwidth A = 
300 ^eV. In the frequency domain, the Lorentzian filter 
has a long tail, Gaussian filter has a shorter tail, and 
rectangular filter has an ideally sharp cut. Therefore, 
from the emission spectrum, the rectangular filter (with 
a bandwidth less than x = 2 meV) may be expected to 
be the most effective filter, but we see in the following 
that the real situation is not so simple. 

In Fig.0](b), we show the (0) obtained for the emis¬ 
sion spectrally filtered by the three types of filters as a 
function of the bandwidth A. As expected, the g^^^O) 
is reduced if the bandwidth is chosen as A < x for all 
filters. On the other hand, if the bandwidth is chosen 
too small, ( 7 ^^^( 0 ) increases as A decreases due to the in¬ 
creased time uncertainty {At = Tc = 1/A in Fig. [3]) as 
mentioned above and in previous literature [2lj. There¬ 
fore, by considering the two opposing effects, spectral 
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emission energy m-(o^ (fieV) fiiter bandwidth i (//eV) XX binding energy / (meV) 

FIG. 4: (a) Emission spectrum of a QD SP emitter (solid) for XX binding energy x = 2000 /reV is shown with the normalized 
filter functions 1^(0;)(dashed) for Lorentzian (black), Gaussian (green), and rectangular (magenta) filters with A = 300 /reV. 
(b, c) g^^\0) of the QD emission after spectral filtering by the three types of filters, Lorentzian (black), Gaussian (green), and 
rectangular (magenta); (b) hlter bandwidth (A) dependence for XX binding energy x ~ 2000 /reV, and (c) The x dependence 
of the g^^\0) at the optimal filter bandwidth Aopt (~ 100 geV for Lorentzian filter in (b); gr*'^^(0) at Aopt are also indicated in 
(b)). We set (xsp, Fph) = (0.67, 20) geV, small pump rate P in the linear regime 0, and the spin flip time ts = xs = 10 ns 
for all figures. 


suppression of the unwanted detection of XX emissions 
and increasing time uncertainty (decreasing time resolu¬ 
tion) for too narrow filters, the existence of the optimal 
hlter bandwidth Aopt is expected. As predicted from the 
above argument, we found the ^^^^(0) shows the minima 
(=0.0027 for Lorentzian, 0.0025 for Gaussian, and 0.0048 
for rectangular hlters) in Fig. 0] (b) (Fph < Aopt < X 
where Fph gives the exciton linewidth). Here, we note 
that our results for Lorentzian hlter perfectly agree with 
those obtained by using the other method |15l |. numeri¬ 
cally showing the validity of our method. 

The above hndings, e.g. the existence of an optimal 
hlter bandwidth at Fph < Aopt < X; seem to be trivial. 
However, following hndings are rather counter-intuitive; 
(I) the rectangular hlter has the largest minimum value 
of (0) among the three hlters although the rectangu¬ 
lar hlter ideally cuts the XX emission in the frequency 
domain, and (H) the optimal hlter bandwidth Aopt is 
much larger than the emission linewidth ^ Fph. (I) also 
applies to a wide range of the XX binding energy (0.5 
meV< X < 8 meV) as seen in Fig. 0] (c), where the min¬ 
ima of g^^\0) as a function of x are shown for the three 
hlters. Fig.0](c) shows that the Gaussian hlter will be the 
best hlter to purify the SP emission from this neutral QD 
system (after the optimization of the bandwidth). (I) can 
be understood as the difference of the hlter correlation 
function in the time domain, /(r). The /(t)(= frij)) in 
Eq. (15^ is the sine function with the slow power-law de¬ 
cay at large r, different from the fast exponential decay 
for the other two hlters. Therefore, the increase in the 
time uncertainty matters signihcantly if a rectangular hl¬ 
ter is used. In the case of the Gaussian hlter, because the 
correlation function is Gaussian also in time domain, the 
long-time tail is strongly suppressed compared with the 
rectangular hlter. Therefore, the lower value of g^'^^ (0) in 
Fig. 0] (c) with the Gaussian hlter is reasonable. 


B. SP emitters under coherent pumping 

Our next example to study the hltering effect is a 
resonantly driven SP emitter. The resonantly scattered 
light by an emitter exhibits Mollow triplet emission spec¬ 
trum [^, which can be applied to an indistinguishable 
SP source [1, [Tl|, [H, [2l|, since the resonant excitation 
prevents emitters from suffering spectral diffusion and 
dark exciton effects, and by reducing dephasing pro¬ 
cesses. 

The physics of the SP emission from the scattered light 
is illustrated in Fig. 0] (a). In the presence of a coherent 
laser held (frequency uig) in resonance with the emitter 
{ujx = wl = wq), the scattered light is known to ex¬ 
hibit the Mollow triplet [2^ (Fig- [5] (b)) with a central 
(G) peak and two side peaks (L, U). The three peaks in 
the huorescence spectrum correspond to the transitions 
(arrows) between the dressed states indicated by C (dash- 
dotted), L (dashed), U (solid), respectively in Fig. [5] (a). 
From the illustration, the side peak (say the upper, U) 
is successively followed by emissions of the other side (L) 
or central (G) peaks. Therefore, successive two-photon 
emission within the same side peak is strongly suppressed 
if the splitting between the dressed states is larger than 
the line width (2r2R ^ if the linewidth is limited by 
the emitter lifetime). This scheme to produce highly ef¬ 
ficient and distinguishable single photons has been stud¬ 
ied in recent years with QD emitters @,[11 [11 mi , in 
which the spectral filtering of a side peak emission (say 
Lo = luq + 2Dr) is essential. In this scheme, the major 
cause of the contamination noise on the SP purity is the 
other two emission peaks (w = wq, ujq — 217^) and the 
excitation laser itself (w = wq). The main physics can be 
described by a resonantly driven two-level system, 

4-P = iipH] +jspCa^p+ {Tpi,/2)CaJ, (39) 










bare states dressed states 



FIG. 5: Illustration of the single photon emission from coherently driven emitters (the Rabi frequency fin, and wo = = ^x, 

laser frequency uil, and the emitter transition frequency uix)- (a) the bare and dressed energy levels with the radiative 
transitions marked by arrows corresponding to the emission peaks in (b), and (b) the fluorescence spectrum showing Mollow 
triplets, a central peak (C) &t lj = loq and upper (U) and lower (L) side peaks. The of the upper side peak emission 

(U) spectrally selected by three types of filters with ljJf = cuo + 2flR, Lorentzian (black), Gaussian (green), and rectangular 
(magenta), (c) 3^^^(0) as a function of the normalized filter bandwidth A/Hr,. (d) g^^\0) minimized in the range 0 < A < 2flR 
as a function of the normalized rabi energy IIr/7sp (the values inside brackets and arrows indicate the minima in (c), the upper 
limit of the range is set in order to avoid the detection of the driving laser light in the shaded area, A > 2f2R). We set the 
spontaneous emission = O.SOr and dephasing rate Fph = 0 for (b) and (c), and the dephasing rate Fph = 0 for all hgures. 


where H = 17^(7“'" + in the rotating frame, and 

'Ysp and Fph are the spontaneous emission and dephas¬ 
ing rates. Here we set the dephasing rate Fph = 0 
according to the experimental reports showing lifetime 
limited linewidth [l2|. For the calculation of the emis¬ 
sion properties, we define the emission field operator by 
= (7^, and the central frequency of the filter wf is 
set ujp = uiq -f- 217^. 

Figure [5] (c), shows the simulated ^^^^0) of the upper 
side peak emissions spectrally filtered by the three types 
of filters, Lorentzian (black), Gaussian (green), rectangu¬ 
lar (magenta) filters. The (0) is plotted as a function 
of the filter bandwidth A, which shows a minimum in the 
regime 0 < A < 217^ reflecting the physics of SP gener¬ 
ation mentioned above; If A < 2f7R, g^^^(O) decreases as 
A decreases since the spectral selection for the side peak 
becomes effective. On the other hand, for A much less 
than the spontaneous emission rate jsp or the linewidth 
( 7 sp = O.SHr in Fig. M (c)), increases as A de¬ 

creases due to the degraded time resolution. The 
at the minima depends on the types of the spectral filters, 
being similar to Fig. 0] (b), the result for an incoherely 
pumped QD. In this case, however, the Gaussian filter 
gives the smallest i-O- the most pure SP emis¬ 

sion, and the Lorentzian filter gives the worst purity for 
'ysp = O.SHr (2 f7R/7sp = 6.7) in Fig. [S] (c). Here, we 
again note that our results for Lorentzian filter perf ectly 
agree with those obtained by using the other method , 
numerically showing the validity of our method. 

In Fig. [S] (d), the g^^^(O) at the minima is shown as 
a function of the ratio between the Rabi frequency over 
the linewidth, 2I7 r/7sp. It is clearly found that g^^^(O) 


at the minima decreases (i.e. the SP purity increases) 
as 2I7 r/7sp increases. This is because the contamination 
source i.e. the other emission peaks, G and L, become 
spectrally separated and suppressed well by the filters for 
larger splitting, 2f7R/7sp. An interesting feature is that 
the Lorentzian filter gives larger than others i.e. 

the performance to obtain high SP purity is the lowest 
among the three types, while the rectangular filter was 
the worst choice in the case of incoherent excitation in 
Fig. |4] (c). Moreover, an interesting result is that the 
most efficient filter type to give the highest SP purity 
depends on the pump parameter 2I7 r/7sp in Fig. [S] (d). 
For the Rabi splitting not too large 2I7 r/7sp < 60, the 
Gaussian filter with the smallest g*'^^(0) is the best choice 
among the three. For the strong Rabi field 2I7 r/7sp > 
60, the rectangular filter with the smallest g^^^(O) is the 
best filter type. We should note here that the latter case 
especially is quite different from the result of incoherent 
excitation in Fig.|4](c). 


A simple explanation for this nontrivial matching be¬ 
tween the filter and emitter remains elusive. In general, 
the matching between filter types and emitters will de¬ 
pend on the non-universal details of the pumping and 
emission dynamics. Therefore, the numerical simulation 
under given conditions and the direct comparison be¬ 
tween different filters are necessary to tailor and optimize 
the quantum emissions. 
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IV. SUMMARY 

We proposed a calculation method, based on the super¬ 
operator eigenvalue decomposition technique, for photon 
statistics of spectrally filtered fields with various types 
of filters. This method can give exact results when the 
emission dynamics is given by quantum master equations, 
which can be applied to a wide variety of quantum emit¬ 
ters, and solvable with the eigenvalue approach (matrix 
diagonalization). Also, it is possible to treat a wide vari¬ 
ety of filter function if analytic expressions for the convo¬ 
lution functions, s and Zk in Eq. (m and Eqs. (|T6|) - (|T8|) . 
are obtained. 

As typical examples, focusing on three filter types, 
Lorentzian, Gaussian, and rectangular filters, we applied 
this method to QD single photon (SP) emitters. With 
the simulation for two cases, under incoherent excita¬ 
tions and under coherent and resonant excitations, we 


found condition-dependent non-universal matching be¬ 
tween filter-types and emitters in order to have the high¬ 
est SP purity. 

An interesting issue remaining will be extending this 
method to simulations for periodic and short-pulsed 
pumping, which will allow us to study the effect of the 
spectral filtering in case of short-pulse excitations aiming 
at more realistic operations of the QD SP source (l^ . 
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